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Entangling and disentangling capacities are the key manifestation of the nonlocal content of a 
quantum operation. A lot of effort has been put recently into investigating (dis)entangling capacities 
of unitary operations, but very little is known about capacities of non-unitary operations. Here we 
investigate (dis) entangling capacities of unital CPTP maps acting on two qubits. 



I. INTRODUCTION 

Entanglement content is one of the fundamental ways 
to characterize nonlocal quantum resources (nonlocal 
states and operations). For pure bipartite states the ul- 
timate measure of entanglement, the von Neumann en- 
tropy of entanglement, had been recently discovered 
A universal measure of entanglement for mixed states 
had not been found yet and different measures are used 
depending on the operational context. Nevertheless, the 
important feature of all entanglement measures of states 
is that their values are directly inferred using the param- 
eters of a state itself. 

Similarly to mixed states, the entanglement content 
of quantum operations can be characterized in different 
ways, e.g. via the amount of entanglement necessary to 
generate that operation or via the amount of entangle- 
ment the operation is able to produce/destroy (the so 
called entangling/disentangling capacities). This article 
is concerned with the two latter measures. 

Unlike the amount of entanglement in a state, the 
(dis)entangling capacities of an operation do not have 
an operational interpretation on their own. They mani- 
fest themselves via the change of the entanglement of a 
particular state that the operation acts upon. And the 
operation has to act on a specific state (the "optimal" 
state) in order to realize its (dis) entangling capacity in 
full. Thus, the straightforward way to calculate these 
quantities is to maximize the change of entanglement over 
all possible initial states. 

Substantial progress have been made recently in inves- 
tigating (dis)entangling capacities of unitary operations. 
The capacities of two-qubit unitary operations were ex- 
plicitly calculated 0, H, H| • It was also shown that single- 
shot capacities are equal to asymptotic capacities [H, U 6 1 . 
Some results for higher dimensions were also obtained 7 1 . 
However, extending these techniques to systems of higher 
dimensionality seems to be a very difficult task. Even in 
the two-qubit case the capacities of a general unitary op- 
eration have been calculated numerically, no analytical 
technique is known. 

In all real situations an experimentalist never deals 
with a perfect unitary in the laboratory. And it is need- 
less to say that calculating capacities of non-unitary op- 
erations, i.e. nonlocal quantum maps, is even a bigger 
challenge. 

In this article we consider nonlocal completely positive 



trace preserving (CPTP) maps of the form 

k 

where Uk are unitary transformations. Maps of this type 
are often called random unitary processes, and they are 
doubly stochastic. The map ([1]) may arise, for example, 
if the desired unitary transformation can be implemented 
successfully only with certain probability, while another 
unitary is realized in the case of failure. A continuous 
version of the map |T]) may arise naturally in experiment 
if parameters of a desired unitary transformation are sub- 
ject to a noise (the case of Gaussian noise will be analyzed 
in detail in Section llV Cp . The scope of this article covers 
the case of r that act on two qubits. We calculate single- 
shot (dis)entangling capacities of r in some particular 
cases. 

The structure of the article is as follows. In Sec. |TT] 
the definition(s) of (dis)cntangling capacities of unitaries 
are presented and some recent results concerning two- 
qubit unitaries. Sec. IIIII generalizes the definition of 
(dis)entangling capacities for non-unitaries. Some nu- 
merical results for (lower bounds on) (dis)entangling ca- 
pacities for discreet and continuous mixtures of unitaries 
are presented in Sec. [TV] 

II. (DIS)ENTANGLING CAPACITY OF A 
UNITARY: DEFINITIONS AND SOME RELATED 
RESULTS 

Consider an unitary operation Uab that acts on a 
tensor product Hilbert space Ha <8 Hb of two spatially 
separated particles A and B. If Uab can not be de- 
composed into a tensor product of local unitaries, i.e. 
Uab 7^ Va® Wb, then we say that Uab is nonlocal. Un- 
like local unitaries, nonlocal unitaries have an ability to 
produce or destroy entanglement. This ability is usually 
characterized by the entangling, E'(U), and the disen- 
tangling, E^(U), capacities, i.e. by the maximal increase 
(decrease) of entanglement that can be achieved when 
U acts on quantum states. To quantify these capacities 
we have to choose appropriate measures of entanglement. 
The most sensible choice is to use the entanglement of for- 
mation [Tll | as a measure of entanglement of the initial 
state p, and the distillable entanglement [l2| as a measure 
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of entanglement of the final state UpU^ . The reason for 
this asymmetric choice is purely operational one. What 
counts is the amount of resources (pure maximally en- 
tangled states) needed to create p (asymptotically) and 
the amount of pure-state entanglement one will be able 
to extract from U pW , again asymptotically. Thus the 
most general definition is 

E\U) = max[D{UpU^) - E F (p)], 

EHU)=max[E F ip)-D(UpU% 
p 

where the maximization is over all possible states p 
(mixed and pure) accessible to U . The Hilbert space of 
an accessible p is not necessarily restricted to Ha ®Hb- 
It turns out to be the case that some U create more en- 
tanglement if the original particles are entangled with 
local ancillary particles 0, y| . It also appears to be the 
case that the maximization in Eq. @ can be restricted 
to pure-states only therefore, the definition ((2|) can 
be simplified as 

E\U) = m&x[E{U\4>)) - E(\ip))} 

(3) 

E l (U) = max[E(\tp)) - E(U\iJj})}, 

ili 

where E is an entanglement measure for pure state 
(Throughout this paper we will use the von Neumann 
entropy of entanglement as the most appropriate mea- 
sure). This obviously simplifies the job significantly. 

Let us briefly recall the main results for A and B being 
two-level particles, qubits. 

Any V ab acting on qubits can be decomposed as (2l.[Toj 

U AB = [V A <8> V B ] e*^™.' [W A ® W B ] , (4) 

where 7r/4 > £ x > £ y > \£ z \ > 0. The middle term sand- 
wiched by local unitaries is called the canonical decom- 
position of U. Any U can be transformed to its canon- 
ical form by sandwiching it with Hermitian conjugates 
of corresponding local unitaries. That means that the 
canonical form is genuinely nonlocal part of U - every- 
thing else is local. The beauty of this results is that 
out of 15 real parameters that parameterize a general 
two-qubit unitary only three are necessary to describe 
its nonlocal nature. It simplifies considerably the clas- 
sification of nonlocal unitaries. For the purpose of our 
discussion three classes can be identified; namely, the 
Controlled-NOT(CNOT)-class ^ 0, & = & = 0), 
the DoubleCNOT-class ^ Q, \ y ^ 0, f, = 0)j_and the 
S WAP-class (all three £ Q are not equal zero) [H, [l6| . The 
names reflect the fact that the corresponding "mother" 
unitary transformation (i.e. with £ a = 7r/4 for a =/= 0) 
belongs to that class. 

The main results for qubits are @, 0] : 

(a) El(U) =E±(U). 

(b) For CNOT-class the optimal state, i.e. the state 
that satisfies definition ([3]), lives solely in the Hilbert 



space of particles A and B (no ancillas are needed) and 
takes the form 

\r pt ) =cosa|0)A|0) B ±isina|l) A |l) B , (5) 

where ± correspond to E' and E- 1 respectively. Thus 
all U from that class achieve their capacity by acting on 
pure states with the same Schmidt basis (only values of 
Schmidt coefficients differ depending on the value of ^ x ). 
The values a = f(£, x ) can be obtained by straightforward 
numerical optimization. 

(c) If £ Q < 7r/4 the maximal capacity is achieved when 
\ip opt ) is already entangled. 

(d) Unitaries of the CNOT-class achieve their capaci- 
ties by acting on optimal states that lie in H.a®T~(-b ■ How- 
ever, unitaries of the DCNOT and SWAP-classes achieve 
their capacities only if the original particles are entangled 
with local ancillas. It was conjectured that it is sufficient 
to take the size of ancillas equal to the size of original 
particles. This conjecture was supported by numerical 
simulations for qubits [3[. 

III. ENTANGLING AND DISENTANGLING 
CAPACITIES OF A NON-UNITARY 

For non-unitaries we will use a definition similar to Eq. 
©, i.e. we define 

E^t) = miac[D(T(fi)) - E F {p)] 

(6) 

E l {T)=mzx[E F (p)-D(T(p))}. 
p 

However, in general here we cannot justify reducing 
the search to pure states. This is due to the fact that 
distillable entanglement is not necessary a convex mea- 
sure. 

We can argue, nevertheless, that in the case of mix- 
tures of unitaries acting on two qubits without ancillas 
the distillable entanglement can be regarded as a convex 
measure. Indeed, a mixture of optimal states ([5]) forms a 
Bell-diagonal state for which the lower and upper bounds 
on distillable entanglement [TtJ 

S{ PA ) - S(pab) < D(p) < E re { P ab) (7) 

coincide. Here we recall that the relative entropy of en- 
tanglement Ere(x) is a convex measure. 

If ancillas are used then the situation is more compli- 
cated. We leave the question of whether the capacities 
are attained on pure states as an open question and cal- 
culate the lower bounds on these capacities using pure 
states. 

IV. MIXTURES OF UNITARIES ACTING ON 
TWO QUBITS 

The properties of two-qubit unitaries described above 
m Section HJ can help us to generalize that approach to 
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FIG. 1: (color online) E^ (r) (solid line) and E^(t) (dashed 
line) as functions of A for different values of £, where p = 0.5. 
The highest curve corresponds to £ = n/4. The lowest curve 
corresponds to £ = 0. Here A is measured in radians. 

mixtures of unitaries as in Eq. (jTJ) [18]. Here we use two 
methods for calculating (r) and E^(t). 

Method I: We make an assumption about the partic- 
ular form of the optimal input state, and subsequently 
find the optimal values of its parameters. 

Method II: We perform a direct numerical optimiza- 
tion without making any a priori assumption about the 
optimal state (except of its purity). 

A. Example I: Discreet CNOT-mixtures 

Consider a mixture of unitaries of the CNOT-class, 

Uk = e*>*°*. (8) 

Here we use Method I. From continuity it follows that 
the optimal state is expected to lie on the 2-dimensional 
manifold of (superpositions of) states of the form ([5]) or 
their convex mixtures. Moreover, in this special case we 
can adopt the argument of [|[ (see Sec. |TTJ) and claim that 
the search can be restricted to pure states only. Thus, 
the state optimal for r is again of the type ((5|) . 

As a simplest case let us consider only two unitaries 
Ui and C/ 2 : 

7-00 pUipu} + (1 - P )u 2P ul (9) 

For convenience let us define A = £| — and denote 
simply by £, so £^ = £ + A. We will fix £ and analyze 
J57T and for various A. For A = the map reduces 
to a unitary (with an appropriate capacity). As smaller 
angle means smaller E'(U), we would expect that if U\ is 
mixed with U2, where A < 0, then the entangling capac- 
ity of the resulting map will decrease relative to E^(U). 
This intuition is consistent with the results presented on 
Fig. [TJ Similarly, we might expect that the entangling 
capacity of the map will increase with A > 0, and that 



FIG. 2: (color online) DCNOT: E\t) (empty diamonds with 
solid line fit) and E^(t) (empty triangles with dashed line fit) 
as functions of A for & = = tt/8. SWAP: E\t) (filled 
diamonds with solid line fit) and E^(r) (filled triangles with 
dashed line fit) as functions of A for = = = tv/8. In 
both cases p = 0.5. Here A is measured in radians. 

this capacity will reach its maximum for maximal A, i.e. 
maximal E^(U2)- However, Fig. [1] shows that this is not 
the case. E*(j) indeed grows while A is positive and 
relatively small, reaching its maximum for certain inter- 
mediate positive value of A and then starting to decrease. 
In other words, if U\ and p are fixed, then maximal E^ (t) 
is achieved for some intermediate U2 with Q. > but 
not for Q. = 7r/4. This result might seem counterintuitive 
from the first sight, but it has a clear explanation. Let a\, 
«2 be the corresponding optimal values of a in Eq. ([5]) 
for U\ and U2 respectively, then the optimal value of a T 
will lie somewhere in between, i.e. satisfy a\ > a T > a 2 . 
For Q = 7r/4, U2 can realize its entangling capacity of 1 
ebit if q;2 = 0. However, when U2 acts on a state with 
a T > a 2 , then it creates less than 1 — H [(cos a T ) 2 ] ebit. 

Disentangling capacity, E^(t), behaves differently. It 
is monotonic with A. It equals E' for A = as expected, 
and it is strictly larger than E* for all other values of A. 
The last observation shows behavior completely opposite 
to that of unitaries. In a sense, it is easier for non-local 
map to destroy entanglement rather than create it, while 
for unitary operations the opposite holds [7[ . 

This approach can be similarly applied to any finite 
number of unitaries and to continuous distribution of uni- 
taries, which will be discussed in Section [TV CI 



B. Example II: Discreet DCNOT and 
SWAP-mixtures 

For mixtures of unitaries of DCNOT and SWAP-class 
Method II was used. The details of numerical calcula- 
tions are presented in Appendix. We conjectured that 
it is sufficient if local ancillas are qubits. Selected re- 
sults are shown in Fig. [21 We can see that for DCNOT- 
mixture the behavior of E^ (r) and E* (r) is qualitatively 



4 



similar to CNOT-mixture (Fig. []}. However, for SWAP- 
mixture slightly different behavior is obtained. In partic- 
ular E^(t) exhibits maximum at an intermediate value 
of A. It is also noticeable that for D > 0.357, E^ (t) is 
smaller for SWAP-mixture than for DCNOT-mixture. It 
is a counterintuitive result that SWAP-mixture which is 
naturally considered as "stronger" that DCNOT-mixture 
has lower entangling capacity. However, again similar to 
the line of thought in Example I we can argue that for 
(relatively) large A the second unitary, J7 2 is too strong, 
and therefore when it acts on the optimal state (optimal 
for the mixture, not for itself) it causes more destruc- 
tion that corresponding f/ 2 of DCNOT-class would have 
caused. 



Example III: Entangling capacity of noisy 
unitary with Gaussian fluctuations 
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FIG. 3: (color online )_ E T (r G ) (solid line) and E l (t g ) (dashed 
line) as a function of £ for several values of Q: 0.01, 0.18, 0.35, 
0.52, 0.69, and 0.86. The dotted line corresponds to fl = 0, i.e. 
a unitary transformation. £ and Q are measured in radians. 



So far we analyzed discrete mixtures of unitaries. In 
this section we analyze a continuous distribution, which is 
usually what experimentalists deal with. These distribu- 
tions arise due to uncertainty in one or more parameters. 
Such uncertainties may be caused by the limits of cali- 
bration precision of the devices and by high sensitivity of 
systems used to generate desired interactions. For exam- 
ple, the strength of exchange coupling between donors in 
silicon based solid-state architectures for quantum com- 
puting exhibit significant uncertainty resulting in error 
in gate operation 

In particular, we consider the case when a non-unitary 
map arises if a unitary from CNOT-class is subject to a 
Gaussian noise. 

Recent work [20] analyzed the capability of noisy 
Hamiltonians to create entanglement. In particular, in- 
teractions of the form Eq. ((HJ, where £ is Gaussian dis- 
tributed with the mean £ = tt/A and standard deviation 
fi, were considered. Without noise this operation (CNOT 
operation) is able to create a maximally entangled state 
if it acts on a disentangled pure state. The authors an- 
alyzed the situation when the noisy operation acts on 
initially disentangled state, which is by itself subject to 
a Gaussian noise. Its capability to create entanglement 
was characterized in terms of the condition for insepara- 
bility of the resulting mixed state (via Peres-Horodecci 
separability criterion) . 

The aim of our analysis is different. We consider noisy 
interactions with £ £ [0, 7r/4] and calculate their entan- 
gling and disentangling capacities in terms of £ and Q. 
Thus, we give a comprehensive quantitative characteriza- 
tion of the non-local content of these noisy maps in terms 
of their entangling and disentangling capacities. Unlike 
[20l | we do not test the resulting state on inseparability, 
rather calculate its distillable entanglement explicitly. 

The action of a unitary U = exp[i^a^a^], where £ is 
Gaussian distributed with the mean £ and the standard 
deviation fi, on the state pab can be seen as a non- 



unitary CPTP map 

r G (p)^-=-/ e-^UpUUt (10) 

which is a continuous mixture of unitaries of the CNOT- 
class. Similarly to the Sec. IIV Al we consider pure initial 
state, i.e. p = \tp)(ip\, where ijj takes the form ([5]), cal- 
culate the distillable entanglement of the output mixed 
state, and maximize it over a. Figure [3] presents nu- 
merical results for E'(t) and E^(t) as functions of £ 
for several values of fl. We can see that already for 
m 0.01 rad we obtain only very small deviation from 
the (dis)entangling capacity of the unitary, i.e. tq with- 
out noise -O = 0. As 51 increases the disentangling 
capacity increases and the entangling capacity decreases. 
The former fact should not be surprizing as it is known 
that entanglement can be destroyed even by local CPTP 
unital maps [2l( . Thus, the more dispersed the distribu- 
tion of £ becomes, the easier for tq to destroy entangle- 
ment and the harder to create it. Nevertheless, we see 
that even when f2 is relatively large tq is still able to 
create considerable amount of entanglement. 

V. DISCUSSION AND CONCLUSION 

We have discussed the entangling and disentangled ca- 
pacities of nonlocal CPTP unital maps, i.e. maps that 
can be represented as probabilistic mixtures of unitaries, 
and have calculated these capacities in some particular 
cases for two qubits. Three classes of unitaries were con- 
sidered, namely the CNOT, DCNOT, and SWAP classes. 

We have observed that the disentangling capacity was 
always larger than corresponding entangling capacity, 
which contrasts with the unitary case where the both 
capacities are equal for qubits and for higher dimensions 
disentangling capacity cannot be greater than entangling 
capacity |7|. 
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In the case of the CNOT-class our results were 
obtained via straightforward generalization of the 
method for CNOT-class unitaries. We argue that the 
(dis)entangling capacity is achieved when a map acts on 
the optimal pure state from the same family as in the uni- 
tary case. Both discrete and continuous mixtures were 
analyzed. In the case of the DCNOT and SWAP-class 
direct numerical optimization was performed. We have 
conjectured that dimensions of the local ancillas are equal 
to the dimensions of the original particles, i.e. the ancil- 
las were taken to be qubits. 

A number of open question can be addressed in a future 
research. 

It will be interesting and useful to prove (or disprove) 
the general conjecture that the sizes of local ancillas can 
be taken equal to the sizes of original particles. 

Here we have calculated single-shot capacities. In the 
case of unitaries it had been shown that in the asymptotic 
regime one cannot do better @, @. It is important to 
check whether this result holds in the non-unitary case. 

In the case of DCNOT and SWAP-mixtures we per- 
formed maximization over pure states only thereby ob- 
taining lover bounds on E' (t) and (r) , but not their 
actual values. 

In our future research we will address the question of 
whether these bounds are tight. It might be the case 
that optimal states for these operations are mixed and, 
consequently, the capacities are higher than we have cal- 
culated. 
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APPENDIX A 

We have used two-dimensional ancilla on each side. 
Consider a general state of four qubits in the tensor- 



product of the computational bases of the original parti- 
cles A, B and the ancillary particles A' , B' 

\^)AA'BB' = c «J,fc,i|«)A|j>A'|fc>s|Os'- ( A1 ) 

There are 16 terms in the above superposition with 16 
complex amplitudes c% j k i, therefore \^) can be parame- 
terized using 30 real numbers (if we take into account the 
global phase and normalization) . We will parameterize it 
in the following way [22|. First, to facilitate our analysis 
it is easier to incorporate four indexes i, j, k and I, each of 
which runs from to 1, into a single index, x, that runs 
from 1 to 16. This can be done by using the formula 
x = 8i + 4j + 2k + I + 1, which is essentially a formula 
for converting a number from the Boolean representation 
to the decimal. Second, we present amplitudes c x in the 
form 

c x = \c x \e ie *, (A2) 

where Yl x =i \ Cx \ 2 = 1 an< ^ ^ = 0- Third, we introduce 
new parameters <p x such that 
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c x = sin0 x _i Y[ cos 4> y , (A3) 

y=x 

where <fio = tt/2. Thus the state |^) is parameterized by 
30 angles. The advantage of this parametrization is that 
we restrict their values only to the interval [0, 2tt] that 
simplifies numerics. 

We proceed as follows. A program generates a vec- 
tor of 30 random numbers in the interval [0, 2ir}. This is 
the initial state. We then apply the non-local map and 
obtain a final state. We calculate the value of the gain 
in entanglement AS = S(t(^) BB ') - S(t(^) A a>bb') - 
S(Ttaa' 1 1 3/)(\E f | ) . After that we vary the values of the ran- 
dom vector by a small amount and repeat these calcula- 
tions again, thereby obtaining a gradient of the change in 
entanglement in that point. We move along the gradient 
to obtain the next |\?), and the procedure is repeated. 
Eventually, the program reaches the maximum where it 
stops. 
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